Citation: Cell Death and Disease (2012) 3, e397; doi:10.1038/cddis.2012.140 

© 2012 Macmilian Publishers Limited All rights reserved 2041-4889/12 



www.nature.com/cddis 

Proliferation and tissue remodeling in cancer: 
the hallmarks revisited 

EK Markert 1 , AJ Levine 1 ' 2 ' 3 and A Vazquez*' 1 ' 2 ' 4 

Although cancers are highly heterogeneous at the genomic level, they can manifest common patterns of gene expression. 
Here, we use gene expression signatures to interrogate two major processes in cancer, proliferation and tissue remodeling. 
We demonstrate that proliferation and remodeling signatures are partially independent and result in four distinctive cancer 
subtypes. Cancers with the proliferation signature are characterized by signatures of p53 and PTEN inactivation and 
concomitant Myc activation. In contrast, remodeling correlates with RAS, HIF-1a and NFkB activation. From the metabolic point 
of view, proliferation is associated with upregulation of glycolysis and serine/glycine metabolism, whereas remodeling is 
characterized by a downregulation of oxidative phosphorylation. Notably, the proliferation signature correlates with 
poor outcome in lung, prostate, breast and brain cancer, whereas remodeling increases mortality rates in colorectal and 
ovarian cancer. 
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The human body is organized into a hierarchy of level 
structures (cells, organs and organism) with their hierarchi- 
cally interrelated functions. In turn, cancer can be organized 
into its level structures as well: the cancer cell, the primary 
tumor and the systemic disease. On each level, cancer 
deregulates the normal functions and activities associated 
with that level's structure. At the lowest level we have the cell, 
and among its main functions controlled cell renewal and 
proliferation. In turn, deregulation of cell proliferation is a 
prerequisite for carcinogenesis. 1 One level up, cells are 
arranged into tissue, which provides structural support and 
organization, nutrient supplies via blood vessels, and inflam- 
mation and wound-healing activities in response to tissue 
damage or nutrient stress. 1 To form a primary tumor, the 
cancer cells should master the remodeling of the cancerous 
and adjacent tissue, attracting blood vessels and transforming 
the surrounding tissue, to promote tumor growth and 
local invasion. Finally at the organism level, cancer affects 
global functions of the organism, evading its immune 
defenses, utilizing the organism's blood circuitry to spread 
and metastasize, and inducing systemic metabolic alteration 
such as cachexia. 1 It is this final systemic state of the disease 
that ultimately results in patient death. 

This line of reasoning suggests that a systematic classi- 
fication of cancer could display the same hierarchical 
organization. Deregulated proliferation has been studied in 
great detail and we now know a great deal about the signaling 
pathway alterations that allow cancer cells to master this 
major phenotype at the cellular level. 2 ' 3 Tissue alterations, 



such as angiogenesis, hypoxia and inflammation, have 
also been significantly investigated and the signaling 
pathway alterations driving these tissue phenotypes 
discovered. 4-6 Systemic alterations such as metastasis and 
cachexia have been investigated as well, 7 ' 8 albeit to a lesser 
extent due to experimental challenges. These studies have, 
over time, been performed in many cancers, thus generating 
the base for a more universal, across-cancers analysis. 
However, most of the research has focused on the particula- 
rities of individual cancer types that result from the 
overwhelming number of genomic alterations that can drive 
these phenotypes. 9-11 

In this paper, we follow a different approach, focusing on the 
commonalities among all cancers rather than on their 
particularities. To this end, we analyze a large collection of 
gene expression data from various tumor sites. We show how 
cell- and tissue-level aspects comprise two main axes of 
cancer development, proliferation (P) and tissue remodeling 
(R). We leave aside the third systemic level of organization 
because of the present lack of appropriate experimental data 
needed for its characterization. We show how tumors across 
multiple tissues can be characterized by their activities of 
proliferation and remodeling. We furthermore correlate these 
key phenotypes with a variety of the most common pathway 
alterations in cancer. In this context, we also include major 
metabolic pathways and transcriptional patterns reflecting 
loss of lineage or differentiation control, manifested as an 
embryonic stem cell (ESC-like phenotype 12-14 ) or a mesench- 
ymal phenotype. 15-17 Finally, we analyze survival data across 
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cancers, to assess the impact of these two fundamental 
processes as a function of site. 

Results 

Gene expression signatures for proliferation and 
remodeling. Proliferation and remodeling require the expres- 
sion of specific pathways. In turn, the consensus expression of 
genes in these pathways can be used as a surrogate to 
quantify the degree of proliferation and remodeling based on 
gene expression profiles. Gene expression signature analysis 
has been shown to accurately reflect major cellular and tissue 
phenotypes and to provide a robust tool for a systematic 
analysis across data sets. In the case of proliferation these 
pathways include well-known processes associated with cell 
cycle progression and cell division: the G-|/S transition, DNA 
replication, telomere organization, chromosome segregation, 
DNA packaging, the G 2 /M transition and cell division. In the 
category of remodeling, we envision processes involving an 
interaction with the microenvironment or a change relative to 
the microenvironment. Here, we include cell junction organi- 
zation, cell adhesion, cell migration, angiogenesis, cytokine 
production, inflammation and response to wounding. These 
are not by any means exhaustive lists, but a representative 
selection of processes that are involved in proliferation and the 
interaction between the cell and its microenvironment. To 
quantify the manifestation of these different processes, we 
define gene signatures based on the gene ontology annota- 
tions or literature reports, enumerating the genes implicated in 
these functions according to our current knowledge (Materials 
and Methods). 

We tested the resulting proliferation and remodeling 
sub-signatures using reported gene expression data sets 
measured in experimental settings where the desired process 
was activated/inactivated. In an experiment where prolifera- 
tion was induced by IGF-I (Figure 1a), the proliferation 
signatures were indeed present in all samples stimulated 
with IGF-I as opposed to the controls, while the tissue 
remodeling signatures were overall not significantly asso- 
ciated with control or case samples. In order to test the 
remodeling signatures, we chose an experimental setting 



where the epithelial mesenchymal transition (EMT) was 
induced by TGF/? stimulus (Figure 1b). We observe a clear 
transition from blue to red in all tissue remodeling signatures 
starting from the 8-h time point. We also note a decrease in the 
proliferation activity going hand in hand with the onset of 
tissue remodeling, suggesting that cells undergoing EMT 
reduce their proliferation activity while making the transition 
into mesenchymal type. We collected all genes in the sub- 
signatures associated with proliferation and the sub-signa- 
tures associated with remodeling, to create the signatures of 
proliferation and remodeling, respectively. The proliferation 
and remodeling signatures capture the consensus expression 
of the respective sub-signatures in the controlled experiments 
(Figure 1). 

We next investigated the correlation between the 
sub-signatures for the specific cellular processes (e.g., 
DNA replication or cell migration) and the signatures for 
proliferation and remodeling across 3931 samples from 
different cancers, including lung, prostate, breast, colorectal, 
brain and ovarian. The sub-signatures of proliferation are 
highly correlated with the proliferation signature, and the same 
was the case with the remodeling signatures (Table 1). 
However, all proliferation sub-signatures display a negative 
correlation with the remodeling signature, and all remodeling 
sub-signatures are negatively correlated with the proliferation 
signature (Table 1). This corroborates the observations in 
controlled experiments (Figures 1a and b) in the context of 
human cancers. Notably, the negative correlations between 
the proliferation and the remodeling groups are low (bordering 
or lacking statistical significance) compared with the very 
strong correlations within each group of signatures. Together 
these data imply that the two meta-pathways proliferation and 
remodeling measure well-defined and clearly distinguished 
phenotypes; the weak-negative correlation suggests that 
these programs are observed in a statistically nearly 
independent fashion, with a subtle trend toward observation 
of opposing regimes. This is confirmed and further illustrated 
below using a classifying strategy. 

Signaling pathways driving proliferation and remodeling. 

The achievement of the proliferation and remodeling 
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Figure 1 Proliferation and remodeling sub-signatures in controlled experiments. Heatmaps showing the proliferation and remodeling signatures in control experiments. 
The red indicates a significant upregulation, black no association, and blue a significant downregulation of the signature indicated on the right, (a) Induction of cell proliferation 
with IGF-I in MCF7 breast cancer cells, based on data from Creighton et a/. 48 (b) Induction of the EMT with TGF-/? in A549 human lung adenocarcinoma cells, based on data 
from Sartor et a/ 49 
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Table 1 Correlation between gene signatures 



Signature Proliferation Remodeling 



Proliferation sub-signatures 

G^S transition " 0.85 -0.06 

DNA replication 0.87 -0.12 

Telomere organization 0.75 -0.14 

DNA packaging 0.80 -0.17 

Chromosome segregation 0.81 -0.14 

G2/M transition 0.80 -0.14 

Cell division 0.87 -0.06 

Remodeling sub-signatures 

Cell junction organization -0.04 0.66 

Cell adhesion -0.15 0.81 

Cell migration -0.05 0.86 

Angiogenesis -0.04 0.77 

Cytokine production -0.02 0.62 

Inflammatory response -0.07 0.67 

Response to wounding -0.05 0.82 

Signaling pathways 

Myc activation 0.68 -0.12 

p53 inactivation 0.55 0.06 

PTEN inactivation 0.82 - 0.1 1 

HIFIalpha targets 0.20 0.48 

RAS activation 0.10 0.53 

NFKB activation 0.03 0.53 

Metabolism 

Glycolysis 0.36 0.07 

Oxidative phosphorylation 0.33 - 0.28 

SOG pathway 0.47 -0.17 

Stem cell 

Embryonic stem cell 0.80 -0.13 

Mesenchyme development -0.11 0.52 



The Pearson correlation coefficient between the signatures on the first column 
and the signatures of proliferation and remodeling 



phenotypes follows from the activation/ inactivation of signaling 
pathways regulating these processes. In the context of cancer 
this occurs owing to both the activation of oncogenes and the 
inactivation of tumor suppressors regulating proliferation and 
remodeling. To this end, we analyze gene signatures 
quantifying the activation of Myc, Ras, NF/cB and HIF-1a and 
inactivation of p53 and PTEN signaling pathways. Although 
this list is certainly not complete, it covers some of the most 
commonly altered pathways in human cancers. 

We investigated the correlation between the signatures for 
activation/inactivation of signaling pathways and the signa- 
tures for proliferation and remodeling across the 3931 tumor 
samples (Table 1). We observe a strong positive correlation 
between Myc signaling activation, p53 signaling inactivation 
and PTEN signaling inactivation and the proliferation signa- 
ture. This observation agrees with the role of Myc as a major 
promoter of proliferation and p53 and PTEN as major negative 
regulators of proliferation. In contrast, the Ras, NFtcB and 
HIF-1 a signaling activation signatures are strongly correlated 
with the remodeling signatures. This observation was also 
expected, given the role of NFkB in inflammation and of 
HIF-1 a in response to hypoxia. Ras activation has been 
commonly associated with increased proliferation. 18 More 
recently, it has been shown that Kras oncogenic mutations 
activate glycosylation, 19 which is known to have a key role 
in invasion and angiogenesis. 20 Our observation does not 



exclude that Ras activation promotes proliferation under 
certain conditions, but rather suggests regulation of remodel- 
ing as its primary role. 

Metabolic pathways fueling proliferation and 
remodeling. Proliferation and remodeling are both fueled by 
cell metabolism. Proliferating cells need energy and precursor 
metabolites to support the biosynthetic processes required to 
produce a new cell. 21 Similarly, remodeling requires energy 
and precursor metabolites to support processes such as cell 
migration and secretion. For example, quiescent fibroblasts 
have high metabolic activity directed to protein synthesis 22 
These requirements may be satisfied in different micro- 
environment conditions affecting nutrient availability, for 
example, hypoxic conditions. To investigate the metabolic 
requirements of proliferation and remodeling across tumors, 
we focus on signatures for energy metabolism. Glycolysis and 
oxidative phosphorylation (OxPhos) are the main recognized 
pathways for energy generation. Recently, we have discov- 
ered a third candidate pathway for ATP production involving 
serine biosynthesis, one carbon metabolism and the glycine 
cleavage system (SOG pathway). 23 Elements of this pathway 
have been found upregulated in breast tumors 24 ' 25 and 
melanoma 26 and it has been shown to be essential for their 
development. Furthermore, most of the genes in this pathway 
are activated by Myc. 23 

We observe that the signatures of glycolysis, OxPhos and 
the SOG pathway are highly correlated with proliferation 
(Table 1), suggesting that tumors with a high proliferation 
signature are metabolically more active. The same conclusion 
was obtained when analyzing a metabolic model accounting 
for the metabolic demands of proliferating cells 23 This model 
predicts that the rate of glycolysis and the SOG pathway 
increase with increasing the proliferation rate. In contrast, 
OxPhos is predicted to first increase but then decrease with 
increasing the proliferation rate 23 The positive correlation 
between the signatures for OxPhos and proliferation can still 
be explained by the fact that the OxPhos signature is 
negatively correlated with the remodeling signature, which is 
in turn negatively correlated with the proliferation signature 
(Table 1). 

Proliferation and remodeling define different stem cell 
phenotypes. The definition of cancer stem cells capable of 
recapitulating tumor development is quite challenging. The 
EMT can generate cells with stem cell properties, such as the 
CD44 + /CD24- antigen phenotype. 15 ' 27 In an independent 
line of research, it has been shown that subsets of breast, 
liver and prostate cancers manifest signatures characteristic 
of ESCs. 14 ' 28 ' 29 To investigate the relationship between 
these stem-cell like phenotypes and proliferation/remodeling, 
we studied gene signatures for ESC and mesenchymal 
development. We observe that the ESC signature is highly 
correlated with the proliferation signatures, whereas the 
mesenchyme development signature is highly correlated 
with the remodeling signatures (Table 1). The association 
between the ESC signature and proliferation has been 
previously noticed. 12 In fact, ESCs undergo very fast 
proliferation during the initial steps of development. 30 ' 31 
The association between mesenchymal development 
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Figure 2 Gene expression signatures by P/R subtype, (a) Gene signatures expression in tumor samples from different cancer types, (b) Overall gene signature 
enrichment across cancer data sets by P/R subtype. In this heat map, a column represents the subset of tumors in the data set indicated in the bottom that were classified as 
the P/R subtype indicated in the top bar. A red color indicates that the signatures listed on the right is significantly enriched among the samples of that particular data set and 
subtype, the brighter the red the more significant, while the blue indicates the absence of a significant enrichment 
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and remodeling also agrees with previous notions. Indeed, 
the remodeling signatures are characteristic features of 
mesenchyme development. However, to our knowledge, 
there have been no previous studies comprehensively 
analyzing the interrelation between ESC features and 
mesenchymal development in cancer. Our analysis indicates 
that these features might be expressions of stem-like 
structures at different tiers of organization. 

Proliferation/remodeling subtypes across cancers of 
different sites. The correlation of the proliferation and 
remodeling signatures across human cancers indicate that 
they represent partially independent processes. To investi- 
gate this relationship in more detail, we investigate the 
occurrence of tumors with different proliferation/remodeling 
phenotypes. To this end, we classify each tumor sample into 
having (P + ) or not having (P - ) a significant upregulation of 
the proliferation signature and into having (R + ) or not having 
(R - ) a significant upregulation of the remodeling signature. 
Putting together the P and R classification, we define 
four possible subtypes: P-/R-, P-/R + , P+/R- and 
P + /R + . In all cancers analyzed (lung, prostate, breast, 
colorectal, brain and ovarian), a significant number of 
samples were found in each of the four subtypes. This is 
illustrated in Figure 2a for the four major cancers by 
incidence (lung, prostate, breast and colorectal) and in 
Figure 2b for all cancers analyzed. 

We observe a significant enrichment of gene signatures 
for Myc activation and p53 and PTEN inactivation in P + (P + / 
R- and P + /R + ) tumors and for RAS and HIF-1 a activation 
in R + tumors (P — /R + and P + /R + ), in agreement with the 
correlation analysis discussed above (Table 1 ). The glycolysis 
signature is highly enriched in the P+ subtypes in most 
cancers. Glycolysis is promoted by both Myc and HIF-1 a. 32 
However, the Myc signature is enriched in the P+ subtypes, 
whereas the HIF-1 a signature is enriched in the R + subtypes. 
Taken together, these data indicate that the upregulation of 
glycolysis correlates more with Myc than HIF-1 a activation. 
The correlation analysis indicated a positive correlation 
between OxPhos and proliferation, but a negative correlation 
between OxPhos and remodeling. Now, looking at the 
different subtypes, we observe that OxPhos is enriched in 
the P + /R- subtype. OxPhos is activated by Myc but 
inhibited by HIF-1 a. Our observations indicate that HIF-1 a is 
the dominant factor regarding the activity of OxPhos, 
in agreement with previous work. 33 The SOG pathway is 
enriched in the P+ subtypes. This is consistent with the fact 
that most of the genes in this pathway are activated by Myc. 23 
We note, however, that the enrichment is stronger in the 
P + /R - than in the P + /R + subtype, suggesting that some 
signaling pathways involved in remodeling may inhibit the 
activity of the SOG pathway. Finally, we recapitulate the 
association between the ESC signature and proliferation 
and the association between mesenchymal development and 
remodeling, across most tumor sites. 

Association with patient survival. To bring the prolifera- 
tion/remodeling picture into a clinical context, we analyzed 
clinical data associated with the different cancers, focusing 
on overall- or disease-specific-survival (Figure 3). The 



current knowledge would suggest that both the P+ and 
R+ phenotypes associate with tumor aggressiveness 
and poor prognosis. A bird's eye inspection of the Kaplan- 
Meier plots in Figure 3 supports this expectation. Indeed, the 
P-/R- is the group with best prognosis either alone or 
shared with the P + /R- or P-/R+ subtypes. There are, 
however, some cancer-specific differences regarding 
whether proliferation or remodeling is the major factor 
determining survival. 

In lung, prostate, breast and brain tumors, the data 
sets analyzed here indicate that proliferation is the major 
prognostic factor. In lung tumors (Figure 3a), there is a 
significant splitting of the survival plots for the P + and P - 
subtypes (P=2.7x 10~ 4 , log-rank test). In prostate tumors 
(Figure 3b), the P+/R+ group has clearly the poorest 
prognosis (P=3.1 x 10~ 5 , P+/R+ versus the rest), followed 
up by the P + /R- (P=1.2x10~ 5 , P+ versus P-). 
In breast tumors (Figures 3c and d), we observe a significant 
splitting of the survival plots for the P+ and P- subtypes 
(P=5.6x10~ 4 and P=1.3x10~ 7 , respectively). In brain 
tumors (Figures 3e and f), the P+ and P- subtypes 
again result in a significant splitting (P=2.9x10~ 9 and 
P=5.2x10~ 3 , respectively). However, in one cohort, the 
survival plot for the P + /R- group overlaps with that for 
the P -/R+ group (Figure 3e), whereas in another one, the 
plots for P + /R+ and P+/R- groups overlap and are 
separated from the P-/R+ plot (Figure 3f). 

The pattern changes in colorectal (Figure 3g) and ovarian 
tumors (Figure 3h), where the data indicate that remodeling is 
the factor determining survival. In both cases, we see a 
significant splitting of the survival plots between the R + and 
R - subtypes (P= 1 .0 x 1 0~ 4 in colorectal and P= 3.7 x 1 0~ 4 
in ovarian tumors, respectively). It has been previously 
reported that reduced proliferation associates with enhanced 
malignancy in colon cancer. 34 Other reports state that EMT 
predicts poor outcome in colon cancer. 35 These associations 
can be explained by our analysis, showing that remodeling 
associates with poor prognosis and proliferation correlates 
negatively with remodeling. To be precise, in colorectal 
cancers without the remodeling signature, increased 
proliferation associates with poor prognosis (P +/R - versus 
P — /R — ). In contrast, in colorectal cancers with the remodel- 
ing signature there is not a significant association between 
proliferation and survival (P + /R+ versus P-/R + ). The 
latter observation underscores the need to simultaneously 
consider both the proliferation and remodeling signatures to 
make any assessment of their impact on survival. 

Discussion 

The analysis and understanding of patterns of deregulation in 
cancer and their correlation with clinical outcome has come a 
long way. Underlying the enormous body of work describing 
individual gene malfunctions, pathway deregulation and 
more general oncogenic or tumor-promoting events, there is 
an organization of recurring themes, which concern the study 
of the cancer cell itself, cancer tissue and tumor formation, 
and systemic disease. The most prominent of these themes 
are proliferation and tissue remodeling. Many individual 
instances of these basic processes have been previously 
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Figure 3 Proliferation-remodeling impact on survival. Kaplan-Meier survival plots stratified by the P/R subtypes in different cancers. The panels (a-h) represent different 
cancer types as indicated on the top of each panel. 



studied and evaluated with respect to their potential to 
induce metastatic disease and cancer death. For instance, 
an association between proliferation and poor prognosis 
has been noticed before in the context of breast cancers, 36 ' 37 
and in renal and lung cancers as well. 38 Other studies 
have reported associations between signatures characteristic 



of remodeling, 
development. 35 



including hypoxia and mesenchymal 



Our comprehensive analysis reveals that proliferation and 
remodeling can indeed be thought of as two distinct meta- 
pathways characterizing cancers across sites; and that these 
programs can have a different impact on patient survival 
depending on the cancer type. In most cancers tested, the 
P + groups have the worst outcome. However, we found that 
in colon and ovarian tissue, the R + groups had a significantly 
worse survival outcome than the P + groups. This might be 
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explained by the high normal prolific activity within the colon 
(crypts) and the ovaries (follicles), as in an already highly 
prolific environment, structural change might have higher 
impact on tumor formation than additional proliferation. 
However, we only point this out as an interesting observation 
that deserves closer inspection. 

In general, our observations suggest that the defining 
hallmarks of cancer as formulated by Hanahan-Weinberg 1,40 
can be organized into the meta-pathways proliferation and 
remodeling. Clearly, sustained proliferative signaling, evasion 
of growth-suppressing signals and replicative immortality are 
hallmarks leading to sustained proliferation. The acquisition of 
pluripotency or ESC phenotype, based on our analysis, goes 
hand in hand with these proliferative hallmarks. On the other 
hand, promotion of inflammation, induction of angiogenesis 
and activation of invasion can all be seen as an effort of the 
tumor to form a structure viable for growth and independent 
of its originally destined organ tissue. This meta-pathway also 
subsumes mesenchymal development or EMT phenotype, 
which is shown here to associate mainly with remodeling. The 
meta-pathways are glued together by two enabling hallmarks: 
'deregulated metabolism' and 'genome instability and muta- 
tion', with which they interact. Both proliferation and remodel- 
ing represent metabolic investment strategies; as such, they 
will of course both require and/or depend on alterations in the 
metabolic regimes providing energy and building materials. 
The ability to satisfy these metabolic requirements in turn 
depends on the tumor structure and microenvironment, 
particularly the level of vascularization. Genome instability 
and mutation, as already pointed out by Hanahan and 
Weinberg, 40 enables many of the symptomatic hallmarks 
associated with both proliferation and remodeling, as was 
documented here by the signatures for Myc, RAS, HIF-1 a and 
NFkB pathway activation and p53 and PTEN pathway 
inactivation. 

Following our own reasoning, we could continue to ask 
whether cancers manifesting the two meta-pathways together 
(P + /R + ) carry some emerging properties that are not 
established by either of the meta-pathways independently. 
These could then shift the disease to the systemic 
level, inducing metastasis, hijacking immune response, and 
transforming the organism metabolism (cachexia). 
Our analysis indicates that the P + /R+ group has the overall 
poorest prognosis. It is tempting to speculate that P + /R + 
cancers represent the most virulent state of the disease, a 
cancer with such strong tissue-independent organization and 
growth momentum, that it can disseminate itself throughout the 
body, inevitably leading to death. It is also possible that at the 
time of death, excluding inflammation and non-cancer-related 
deaths, all patients carry tumors in this super-cancer state. 

Materials and Methods 

Gene signature analysis. The following signatures were constructed 
based on gene ontology terms: GVS transition (GO:000082), DNA replication 
(GO:0006260), telomere organization (GO:0032200), DNA packaging 
(GO:0006323), chromosome segregation (GO:0007059), G 2 /M transition (GO: 
0000086), cell division (GO:0051301), cell junction organization (GO:0034330), 
cell adhesion (GO:0007155), cell migration (GO:0016477), angiogenesis 
(GO:0001525), cytokine production (GO:0001816), inflammatory response 
(GO:0006954), response to wounding (GO:0009611), glycolysis (GO:0006096), 
OxPhos (GO:0006119) and mesenchyme development (GO:0060485). The SOG 



pathway signature was constructed using the gene list reported in Vazquez etal. 
The remaining signatures were obtained from literature reports: p53 inactivation, 41 
PTEN inactivation, 42 Myc targets, 43 ESC 28 ' 44 and HIF-1 a targets, 45 and NFkB 
targets 46 The genes contained in each signature are reported in the 
Supplementary Table 1. For each signature, we used gene set enrichment 
analysis 47 to determine concordant up- or downregulation of all genes in the 
signature compared with the overall mean expression of genes. We use as 
readout heatmaps, where a red color indicates significant positive association with 
the signature ('ON', P+ <0.05), a blue color indicates negative association with 
the signature ('OFF', P-<0.05) and a black color indicates no significant 
association. 

Cancer data sets. The cancer data sets analyzed here are described in 
Supplementary Table 2. 

Clustering. Cancer samples were clustered based on their proliferation and 
remodeling signature. A sample was defined P+ whenever it manifested a 
significant expression of the proliferation signature, P - otherwise. A sample was 
defined R + whenever it manifested a significant expression of the remodeling 
meta-signature, R - otherwise. The samples were then clustered into the four 
possible subtypes P-/R-, P-/R + , P + /R- and P + /R + . 
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